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Abstract 

Two relativistic approaches are considered to evaluate the quasielastic double-differential and 
integrated neutrino-nucleus cross sections. One, based on the relativistic impulse approximation, 
relies on the microscopic description of nuclear dynamics using relativistic mean field theory, and 
incorporates a description of the final-state interactions. The second is based on the superscaling 
behavior exhibited by electron scattering data and its applicability, due to the universal character 
of the scaling function, to the analysis of neutrino scattering reactions. The role played by the 
vector meson-exchange currents in the two-particle two-hole sector is also incorporated and the 
results obtained are compared with the recent data for neutrinos measured by the MiniBooNE 
Collaboration. 

PACS numbers: 25.30.Pt, 13.15.+g, 24.10.Jv 
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I. INTRODUCTION 



The data on muon neutrino charged-current quasielastic (CCQE) cross sections recently 
obtained by the MiniBooNE collaboration and its comparison with several theoretical 
calculations, have led to an important debate concerning the role played by various ingre- 
dients entering in the description of the reaction: nuclear dynamics (final-state interactions 
(FSI), low-lying nuclear excitations, effects beyond the impulse approximation (lA), etc.), 
as well as possible modifications of the single-nucleon form factors. Although no definitive 
conclusions are yet in hand, a detailed study of modeling versus experiment for inclusive 
quasielastic electron scattering and its extension to neutrino processes can shed light on the 
different interpretations of the discrepancy between theory and experiment. 

When a dipole shape is assumed for the axial form factor, the nucleon axial mass Ma can 
be considered to be the only free parameter within the Relativistic Fermi Gas (RFG) model, 
presently used in many Monte Carlo codes employed in the analysis of the experimental 
data. When compared with MiniBoone CCQE data, the RFG underestimates the total 
cross section unless an axial mass of the order of 1.35 GeV/c^ is employed in the 
dipole prescription for the form factor. This value of the axial mass is considerably larger 
than the accepted world average value [2I, thus yielding a larger axial form factor. This 
should be taken more as an indication of incompleteness of the theoretical description of the 
MiniBooNE data based upon the RFG, rather than as a true indication for a larger axial 
mass. 

For instance, although the RFG incorporates a fully relativistic treatment, required by 
the kinematics of the experiment (mean neutrino energy flux, (E^) = 788 MeV, with values 
up to 3 GeV), its description of the nuclear dynamics is clearly too crude to draw specific 
conclusions on the value of the anomalous axial mass from the departure of the RFG from 
experiment, but rather as a hint on the importance of nuclear effects in describing these 
experimental data. 

However, at the level of the impulse approximation, a number of much more sophisticated 
descriptions of the nuclear dynamics other than the one represented by the RFG, based for 
instance on the use of realistic spectral functions 3|-|5j , when compared with the MiniBooNE 
experimental data also underpredict the measured CCQE cross section, in this respect not 
doing a better job than the RFG. One important consideration that must be taken into 
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account is the fact that, even when these models provide a much more reahstic description 
of the nuclear dynamics than the RFG, they are built on non-relativistic approaches that 
are likely questionable at the kinematics of the MiniBooNE experiment. 

Among the difficulties that one faces when comparing models, is that the effect of the 
ingredients in the model, such as interactions in the final state, may differ greatly from 
model to model. For example in j4| the FSI are barely seen, causing only a simple ~10 MeV 
shift of the QE peak. However, relativistic and semi-relativistic models of inclusive QE 
(e, e') reactions which included a relativistic mean field, that is, described the FSI by means 
of strong relativistic potentials or their semi-relativistic equivalents, have clearly shown the 
essential role played by FSI in order to describe properly the behavior of data |6l-ll0|. 

In addition to the relativistic treatment of the nuclear excitations, in some regions of the 
wide range of neutrino energies where the neutrino fiux for the experiment has significant 
strength, the reaction may have sizable contributions from effects beyond the lA. For in- 



stance. 



in 



n 



Il2| when the theoretical results incorporated multiple knockout excitations, 
they were shown to be in accordance with the total cross section data without the need to 
increase the value of Ma- However, no comparison with the experimental double-differential 



cross section is shown in 
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12| . Moreover, these calculations are based on non-relativistic 



reductions whose reliability at MiniBooNE kinematics may be doubtful. In fact, the kine- 
matics of the MiniBooNE experiment demands relativity as an essential ingredient; not 



only relativistic kinematics should be considered, but also the nuclear 



rent operators should be described within a relativistic framework 13|, iM]. Furthermore, 



dynamics and cur- 



the wide range of neutrino energies, at least for some specific conditions, may also require 
one to account for effects not included in models devised for quasi-free scattering. This is, 
for instance, the situation at the most forward scattering angles where a very significant 
contribution in the cross section may come from very low-lying excitations in nuclei 15 |. 

A systematic analysis of the world inclusive (e, e') data has clearly demonstrated that, for 
sufficiently large momentum transfers, at energy transfers below the QE peak the property 
of superscaling works rather well 16l-ll9|. that is, the reduced cross section, when represented 



versus the scaling variable 



20| , is largely independent of the momentum transfer (first-kind 



scaling) and of the nuclear target (second- kind scaling). Moreover, from the longitudinal 
response a phenomenological scaling function has been extracted that shows a clear asym- 
metry with respect to the QEP with a long tail extended to positive values of the scaling 
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variable (larger energy transfers). Assuming the scaling function to be universal, i.e., valid 



for electromagnetic and weak interactions, in 



21 



22I ] CCQE neutrino-nucleus cross sections 



were evaluated by using the scaling function extracted from (e, e') data and multiplying it by 
the corresponding elementary weak cross section. This approach, denoted simply as "SuSA" , 
provides nuclear-model-independent neutrino-nucleus cross sections, but its reliability rests 
on a basic assumption: the scaling function (extracted from longitudinal (e, e') data) is ap- 
propriate for all of the various weak responses involved in neutrino scattering (charge-charge, 
charge- longitudinal, longitudinal- longitudinal, transverse and axial), and is independent of 
the vector or axial nature of the nuclear current entering the hadronic tensor. In particular, 
the SuSA approach assumes the electromagnetic longitudinal (L) and transverse (T) scaling 
functions to be equal. This property, known as scaling of the zeroth kind, is fulfilled by the 
RFG (by construction) and by most models based on non-relativistic descriptions that, in a 
way, factorize the elementary lepton-nucleus amplitude into a lepton-nucleon part and a part 
containing the nuclear effects 23|, l2J]. Within SuSA, this factorization in the elementary 
amplitude propagates even to the cross section, which is then proportional to the elementary 
lepton-nucleon cross section and to the nuclear response, the latter in this approach being 
a universal function. 

However, from the analysis of the existing L/T separated data, after removing inelastic 
contributions and two-particle-emission effects one finds that the "purely nucleonic" trans- 
verse scaling function is significantly larger than the longitudinal one 25[. This has to be 
attributed to a breakdown of the elementary factorization mentioned before, so that the 
elementary lepton-nucleon vertex inside the nucleus is no longer accurately described by the 
one for free nucleons. One must resort to models such as the relativistic mean field ap- 
proach, denoted as RMF, where the relativistic dynamics introduces significant deviations 
of the behaviour of the elementary lepton-nucleon vertex in the presence of strong scalar 
and vector potentials {23]. This breakdown of zeroth-kind scaling present in the RMF seems 
to be favored by the comparison with data 25 |. 



In a recent paper 



15| SuSA predictions have been compared with the MiniBooNE data 



for the double-differential neutrino cross section showing a systematic discrepancy between 
theory and experiment. Inclusion of 2p-2h Meson Exchange Current (MEC) contributions 
yields larger cross sections and accordingly better agreement with the data. However, theory 
still lies below the data at larger angles where the cross sections are smaller. Before drawing 
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definitive conclusions on the anomalous axial mass, it is important to explore alternative 
approaches that have been shown to be successful in describing inclusive QE (e, e') pro- 
cesses. As just mentioned, this is the case for the RMF, where a fully relativistic description 
(kinematics and dynamics) of the process is incorporated, and FSI are taken into account 
by using the same relativistic scalar and vector energy-independent potentials considered in 
the description of the initial bound states. The RMF model applied to inclusive QE (e, e') 
processes has been shown to describe scaling behaviour, and more importantly, it gives rise 
to a superscaling function with a significant asymmetry, namely, in complete accord with 
data 0, 0] • Moreover, contrary to SuS A, where scaling of the zeroth kind is assumed, the 
RMF model provides longitudinal and transverse scaling functions which differ by typically 
20%, the T one being larger. 

The RMF approach has been applied to the description of CCQE neutrino-nucleus cross 
sections [g], 27H29I and it has been investigated with respect to how scaling emerges from 
neutrino reactions, and how the "theoretical" neutrino scaling functions compare with the 



corresponding ones evaluated for electrons (L and T responses) and with the data 28l. l29|. 

Even at the level of the impulse approximation, the zeroth-kind scaling violation in- 
troduced by the RMF approach, as well as the different isospin character shown by the 
electromagnetic and weak nucleon form factors, can lead to significant discrepancies be- 
tween the results provided by SuSA and RMF approaches. Furthermore, effects beyond 
the lA give rise to additional scaling violations in the transverse responses. Thus, a proper 
relativistic description of these effects is needed in order to compare with the data taken by 
the MiniBooNE collaboration. 

The paper is organized as follows: after this introductory section, in the one following 
(Sect, [ni) we present an analysis of the results obtained using three models for the CCQE 
cross sections, while in Sect. IIIII we end by making a few concluding remarks. 

II. ANALYSIS OF RESULTS 

In this section we discuss the results obtained with the different approaches considered 
and compare with the experimental data. Details on the various approaches considered have 
been presented in previous works. In particular, the SuSA approach and its extension to 



CC neutrino reactions can be reviewed in 22], whereas the basic ingredients entering in 



5 



;he RMF model applied to inclusive electron and CCQE neutrino reactions are given in 
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28|-|30| 



We first show results for the CCQE z/„- C double-differential cross section averaged over 



the neutrino flux ^{E^), namely 



a 



dT^d cos 9 



tot 



d^a 



dTf^d cos 9 



^E,)dE,, 



(1) 



where and 9 are the kinetic energy and scattering angle of the outgoing muon, E^, is the 
neutrino energy and ^tot is the total flux. For each value of the neutrino ener gy the above 
cross section can be expressed in terms of seven nuclear response functions as 
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" d^a ' 




dT^d cos 9 


= o-Q 

E, 



VlRl^ + VccRcc + '^^clRcl + ^llRll 
+Vt {RV + ^T^) + Wt'RV'^ 



(2) 



where Vi are kinematic factors and the indices L, C, T, T', V, A refer to longitudinal, charge, 
transverse, transverse- axial, vector and axial- vector components of the nuclear current, re- 
spectively. 

In particular, in the SuSA approach each response function can be cast as 



Ri{q,u) 



qkp 



Rr{q,oj)fi^) 



(3) 



where q qnd u are the transferred momentum and energy, respectively, m^v is the nucleon 
mass, kp is the Fermi momentum, i?*"" are the single-nucleon responses, ip{q,uj) is the RFG 
scaling variable (see, e.g., 20| for its definition) and f^ip) is the so-called superscaling func- 
tion, containing the dependence on the nuclear model. In the SuSA model it is given by a 
fit to the experimental longitudinal (e, e') reduced response function |31 |. 

In the RMF case, the weak response functions are given by taking the appropriate com- 
ponents of the weak hadronic tensor constructed from the single-nucleon current: 



dr(l)Av)J\ 



(4) 



where (pB and (pp are relativistic bound-state and scattering wave functions, respectively, 
and Jy^ is the relativistic one-body current operator modeling the coupling between the 



virtual W and a nucleon 



32| . The bound nucleon states are described as self-consistent 



Dirac-Hartree solutions, derived within an RMF approach by using a Lagrangian containing 
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cr, u and p mesons j33l435l|. The outgoing nucleon wave function is computed by using the 



same relativistic mean field employed in the initial state. This incorporates the FSI between 
the ejected nucleon (proton) and the residual nucleus. 

Finally, concerning the description of MEC contributions, we use the fully relativistic 
model of 
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38 



to 



l36l . |37| . In particular in the 2p-2h sector we use the scheme applied in 
electron scattering, where all many-body diagrams containing two pionic lines were taken 
into account. However, it is important to point out that, within the present approach, only 
the pure vector transverse response, -R^^, is affected by MEC. Effects in the axial- vector 
transverse response, as well as the contribution of the correlation diagrams (performed re- 
cently for (e, e') in 39||), should be incorporated into the analysis before definitive conclusions 
on the comparison with data can be drawn. Work along this line is presently under way. 
MEC contributions to v-^'^C reactions have been computed within a somewhat different 



approach both for charged and neutral currents, in 40|, |4l|, where the effect of MEC in the 
cross section was found to be less than 10%. 

In Fig. 1 we show the double-differential cross section averaged over the neutrino energy 
flux as a function of the muon kinetic energy T^. In each panel the results have been averaged 
over the corresponding angular bin of cos 9. In all cases we use the standard value of the 
nucleon axial mass, i.e., Ma = 1-03 GeV/c^. We compare the theoretical results evaluated 
using the three approaches, SuSA (green line), SuSA+MEC (blue) and RMF (red), with the 
MiniBooNE data [l|. The case of the most forward angles, 0.9 < cos^ < 1, has not been 
considered since, as shown in [15|, models based on quasi-free scattering cannot describe 
properly this kinematic situation where roughly 1/2 of the total cross section arises from 
excitation energies below ~50 MeV. 

The analysis of the results corresponding to SuSA and SuSA+MEC approaches and their 
comparison with data were already presented and discussed at length in 15(|. We showed 
that the 2p-2h MEC increase the cross section, yielding results that are closer to experiment, 
speciflcally, for data up to cos 6* ~ 0.6. At larger angles, the discrepancy with experiment 
becomes larger while, on the other hand, the role of MEC is seen to be less signiflcant, that 
is, the difference between SuSA and SuSA+MEC becomes smaller as the scattering angle 
increases. 

Cross sections evaluated with the RMF model also yield reasonable agreement with data 
for smaller angles, the discrepancy becoming larger as 9 increases. However, some differences 
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FIG. 1: (color online) Flux-integrated double-differential cross section per target nucleon for the 
i/^ CCQE process on evaluated in the SuSA (green line), SuSA+MEC (blue) and RMF (red) 
models and displayed versus the muon kinetic energy for various bins of cos 6. The data are 



from MiniBooNE 



l|. The uncertainties do not include the overall normalization error 5A^=10.7%. 



emerge from the comparison between the RMF and SuSA predictions. As observed, RMF 
cross sections are in general larger than the SuSA ones. In particular, in the region close to 
the peak in the cross section, the RMF result becomes larger than the one obtained with 
SuSA+MEC. This holds especially for large scattering angles. On the contrary, SuSA and 
SuSA+MEC get more strength in the region of high muon kinetic energies. This can be 
attributed to the breakdown of zeroth-kind scaling in the RMF, in contrast to the other 
approaches where it is assumed to be satisfied. An approach based on RMF, but invoking 
zeroth-kind scaling, yields results that are much more similar to the SuSA ones. 

To make such a statement more transparent, we compare the double-differential cross 
sections evaluated with the three models, but for fixed values of the scattering angle and 
muon kinetic energy. The results are presented against the neutrino energy. We have selected 
as a representative situation the case cos 6* = 0.45 (panel in the middle of Fig. 1) and two 
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(GeV) £„ (GeV) 

FIG. 2: (color online) Double-differential cross section calculated for fixed values of the muon 
kinetic energy and scattering angle and displayed versus the neutrino energy. Results presented 
for the three models: SuSA (green), SuSA+MEC (blue) and RMF (red). 

values of T^: 0.35 GeV that corresponds to the maximum in the neutrino-fiux-averaged cross 
section, and = 0.65 GeV, located in the tail. Results are presented in Fig. 2. As shown, 
for = 0.35 GeV (left panel) the three models produce roughly the same response in the 
maximum at E,y ^ 0.6 GeV. However, the strength in the tail for higher neutrino energies is 
much more significant for RMF, being reduced for SuSA+MEC and much smaller for SuSA. 
This explains why the RMF neutrino-fiux-averaged cross section is significantly higher at 
= 0.35 GeV (see Fig. 1). 

The situation is clearly different for = 0.65 GeV (right panel in Fig. 2). Here, SuSA 
and SuSA+MEC cross sections are larger (compared with RMF) even in the region where the 
cross section reaches its maximum. On the contrary, for larger (located in the tail), RMF 
becomes higher. However, notice that for these kinematics the neutrino energies involved are 
much larger than in the previous case, namely, > 0.8 — 0.9 GeV; this corresponds to the 
tail in the experimental neutrino flux whose average neutrino energy is 788 MeV. Hence, the 
main contribution in the averaged cross section comes from the region with smaller values 
of El, where the difference between SuSA (and SuSA+MEC) and RMF is larger. As already 
mentioned, if one does an RMF calculation that respects zeroth-kind scaling, the results 
would be essentially in agreement with those of SuSA. 

For completeness in Fig. 3 we show the fiux-integrated cross section averaged over the bin 
0.4 < cos 6* < 0.5 evaluated within the framework of the relativistic impulse approximation 
(RIA), but with different descriptions for the FSI. We have considered the relativistic plane 
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FIG. 3: (color online) As for Fig. 1 for the bin 0.4 < cos 9 < 0.5, but now showing the results 
evaluated with RPWIA (green), rROP (blue) and RMF (red). 
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FIG. 4: (color online) Double-differential CCQE cross section for ^^C integrated over neutrino 
flux versus the outgoing muon scattering angle for various bins of the muon kinetic energy T^. 
Results are given for RMF (red hues), SuSA (green) and SuSA+MEC (blue). 



wave impulse approximation (RPWIA), that is, switching off FSI in the RMF calculation, 
and the use of the real part of the relativistic energy-dependent optical potential, denoted 
as rROP. As already shown in previous works [3, Q]' these two approaches fulfill scaling, 
but give rise to scaling functions that lack the asymmetry shown by data. Moreover, scaling 
of the zeroth kind is also highly respected because of the minor role played in this case by 
relativistic dynamics in the final state. 
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Results in Fig. 3 show that the RPWIA and rROP approaches are very similar for all 
T^-values, being also in accordance with RMF, although here the maximum in the cross 
section is slightly reduced while strength is shifted to larger values of the muon kinetic 
energy. This is a consequence of the differences introduced in the scaling functions by the 
particular description of FSI and its impact on the relativistic nuclear dynamics and the 
isospin (third- kind) and zeroth-kind scaling violations 

In Fig. 4 we plot the neutrino-fiux-averaged cross section versus the scattering angle at 
fixed (averaged over each bin). For low muon momenta the three models tend to un- 
derestimate the data, improving the agreement as increases. As observed, when added 
to the SuSA results, the 2p-2h MEC yield an enhancement of the cross section whose mag- 
nitude increases for more forward scattering angles. This result holds for each bin in T^. 
With respect to comparison with data, some general comments already made for the SuSA 
and SuSA+MEC results also apply to RMF: the last also underestimates the data at 
large muon scattering angles, particularly for small T^. However, some important differ- 
ences between RMF and SuSA-based models also emerge. Let us comment on the general 
trend followed by the RMF results as functions of cos^ that clearly differ from SuSA and 
SuSA-|-MEC. In the six panels presented in Fig. 4 RMF cross sections are the lowest for 
the smallest values of cos^. As we move to more positive cos^, the RMF cross section 
grows faster, lying above the results corresponding to SuSA+MEC in the intermediate re- 
gion. Finally, for smaller values of the scattering angle, namely cos 6 approaching 0.9, while 
RMF inverts its behavior and decreases very rapidly, SuSA and SuSA+MEC approaches 
to cos 6* = 0.9 show a much softer slope. In fact, this is the region where the discrepancy 
between RMF and SuSA-based models can be better appreciated. It is very illustrative to 
point out that the general shape presented by the RMF cross section as a function of cos 6 
fits perfectly well the shape shown by data, although RMF predictions fall below the data for 
small muon momenta. The different behaviour of the models is partly due to the fact that 
the RMF is better describing the low-energy excitation region whereas, as already pointed 
out, the SuSA model has no predictive power at very low angles, where the cross section is 
dominated by low excitation energies and the superscaling ideas are not supposed to apply. 
For this reason we do not show results for the highest cos 6 bin. 

In Fig. 5 we present the results obtained by integrating the fiux-averaged double- 
differential cross sections over cos 6' (upper panel) and (bottom panel), respectively. In 
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FIG. 5: (color online) Results obtained with SuSA, SuSA+MEC, RMF and RFG models. Upper 
panel: Flux-averaged integrated cross section displayed versus the muon kinetic energy. Bottom 
panel: As for the upper one, but now for the flux-averaged muon angular distribution. 



addition to the three models considered in previous graphs, here we also include for ref- 
erence the predictions given by the RFG. It is interesting to remark that, in spite of the 
clear differences shown by the RMF and SuSA predictions for the double-differential cross 
sections (Figs. 1 and 4), the integrated results almost coincide. On the contrary, the 2p-2h 
MEG effects produce a visible enhancement in the cross section that is closer to the exper- 
imental data. The RFG results lie somewhere between the SuSA/RMF and SuSA+MEG 
predictions. 
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FIG. 6: (color online) Total CCQE cross section per neutron versus the neutrino energy. The 
curves corresponding to different nuclear models (see text) are compared with the flux unfolded 
MiniBooNE data 

To conclude this section, in Fig. 6 we display the total QE cross section per neutron 
obtained in the models discussed above as a function of the neutrino energy and compared 
with the experimental data. Note that here the integration is performed over all muon 
scattering angles (—1 < cos6' < 1) and energies < < E^). 

As observed in Fig. 5, the discrepancies between the various models tend to be washed 
out by the integration, yielding very similar results for the models that include FSI (SuSA, 
RMF and rROP), all of them giving a lower total cross section than the models without FSI 
(RFG and RPWIA). On the other hand the SuSA+MEC curve, while being closer to the 
data at high neutrino energies, has a somewhat different shape with respect to the other 
models, in qualitative agreement with the relativistic calculation of 42| . It should be noted, 
however, that the result is affected by an uncertainties of about 5% associated with the 
treatment of the 2p-2h contribution at low momentum transfers. 



III. CONCLUDING REMARKS 



Summarizing, in this paper we extend the previous work presented in 15| where the 
focus was placed on the use of the phenomenological SuSA model and its extension incor- 
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porating the role played by 2p-2h MEC contributions. Here, our main interest resides in 
the predictions provided by the RMF approach. This model has been successfully applied 
to inclusive QE (e, e') processes where it has been shown to be capable of reproducing the 
specific asymmetric shape shown by the experimental scaling function. On the other hand, 
and unlike the SuSA approach which assumes scaling of the zeroth kind, i.e., equal longitu- 
dinal and transverse scaling functions, the RMF provides a transverse scaling function that 
exceeds by about ~20% the longitudinal one. This result seems to be in accordance with 
the most recent analyses of the L/T separated (e, e') data. Thus, this violation of scaling 
of the zeroth kind has visible effects when proceeding to studies of CCQE cross sections. 
Furthermore, the different isospin content, (namely, violations of third-kind scaling) of the 



electromagnetic and CC weak nucleon form factors should also be carefully considered [26 1. 

In this work we apply the RMF model to CCQE neutrino reactions on ^^C corresponding 
to the kinematics of the MiniBooNE experiment. Results for the fiux-averaged double- 
differential cross sections are compared with data and the predictions given by SuSA and 
SuSA-|-MEC models. Generally speaking, the RMF model underestimates the data espe- 
cially at large muon scattering angles and low muon energies. This was already observed 
with SuSA and to a somewhat lesser extent with SuSA+MEC. However, the specific be- 
havior shown by RMF clearly differs from that of SuSA and SuSA+MEC; the maximum in 
(Pa / d cos 6 dT^ as a function of for various bins of cos^ gets higher for RMF, whereas 
the tail at high is more pronounced for the SuSA-based models. Also, the general trend 
shown by the curve corresponding to the double-differential cross section as a function of the 
scattering angle for bins of T^, clearly differs for RMF and SuSA (SuSA+MEC) approaches. 
Here, it is very interesting to point out that the specific shape followed by RMF predictions 
fits perfectly well the slope shown by data. 

The single-differential cross sections shown in Fig. 5 where the three calculations yield 
very similar predictions, with almost the same shape and underpredicting the data, also 
shows that it is very useful to compare double-differential cross sections as in figure 4, where 
differences among models may be more easily seen. 

To conclude, let us note that, in spite of the discrepancies introduced by the models in the 
double-differential cross sections, RMF and SuSA approaches provide almost identical results 
for the single-differential cross section, this being found to lie below the data. Although 
the inclusion of 2p-2h MEC contributions increases the differential cross section without a 



14 



significant change of tfie sliape, and tfius seems to improve tlie agreement with the data 
as shown in Figs. 2 and 5, it is also seen (in Fig. 2, but more clearly in Fig. 4) that 
the shape of the cross section is best reproduced by the RMF and does not improve with 
the inclusion of the 2p-2h MEC contributions in SuSA. It is tempting to hypothesize that 
addition of the 2p-2h MEC effects to the RMF results would lead to reasonable agreement 
in both magnitude and shape with the experimental double-differential cross section. 

Finally, as shown in Fig. 6, the impact of the 2p-2h contribution on the total cross 
section increases with the neutrino energy, suggesting that the data can be explained without 
the need for a large nucleon axial mass. However more refined calculations taking care of 
correlation currents, MEC effects in the axial- vector channel, etc., should be performed 
before definitive conclusions can be drawn. 
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